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Abstract: The collective field formulation of the Calogero model supports periodic 
density waves. An important set of such density waves is a two-parameter family of 
BPS solutions of the equations of motion of the collective field theory. One of these 
parameters is essentially the average particle density, which determines the period, 
while the other parameter determines the amplitude. These BPS solutions are some- 
times referred to as "small amplitude waves" since they undulate around their mean 
density, but never vanish. We present complete analysis of quadratic fluctuations 
around these BPS solutions. The corresponding fluctuation hamiltonian (i.e., the 
stability operator) is diagonalized in terms of bosonic creation and annihilation op- 
erators which correspond to the complete orthogonal set of Bloch-Floquet eigenstates 
of a related periodic Schrodinger hamiltonian, which we derive explicitly. Remark- 
ably, the fluctuation spectrum is independent of the parameter which determines 
the density wave's amplitude. As a consequence, the sum over zero-point energies 
of the field-theoretic fluctuation hamiltonian, and its ensuing normal-ordering and 
regularization, are the same as in the case of fluctuations around constant density 
background, namely, the ground state. Thus, quadratic fluctuations do not shift the 
energy density tied with the BPS-density waves studied here, compared to its ground 
state value. Finally, we also make some brief remarks concerning fluctuations around 
non-BPS density waves. 
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1. Introduction 

In the limit of infinite number of particles and finite density, the Calogero-Sutherland 
models |l[, [^ may be reformulated in terms of hydrodynamic variables, namely, the 
particle density and current. The hamiltonian of this collective field reformulation is 
given by g] 

p X y ^^^^^ 

where H^ing denotes a singular contribution 



X f P 

^ 2m J ' x-y 



A-1 



dxd,^5{x-y)l,, (1.2) 
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and P is the principal part symbol. 
Here, 



N 



p(x) = 5(x — Xj) (1.3) 



1=1 



is the collective - or density - field, and 



6p{x) 



is its canonically conjugate momentum. It follows from (|1.3|) that the collective field 
is a positive operator 

p{x)>0, (1.5) 
and that it obeys the normalization condition 

dxp{x)=N. (1.6) 



The latter constraint is implemented by adding to (|l.lj ) a term pi J dx p{x) — N 

where p is a Lagrange multiplier (the chemical potential). 

The first term in (|1.2| ) is linear in p{x). Therefore, its singular coefficient 

— ^dx amounts to a shift of the chemical potential u by an infinite con- 

y=x 

stant. The last term in ( |1.2|) is, of course, a field independent constant - an infinite 
shift of energy. 

For the sake of being self-contained, we briefiy explain in the Appendix how to 
derive the collective hamiltonian ( |1.1[ ) and its singular part ( |1.2| ) from the microscopic 
Calogero hamiltonian. 

It is worth mentioning at this point that the Calogero model enjoys a strong- 
weak-coupling duality symmetry 0. At the level of the collective Hamiltonian 
( |lTD , these duality transformations read 



A ' " ~X ' " ~ A ' '^^^^ " -^PiV ) and n{x) = — , (1.7) 

and it is straightforward to see that these transformations leave (|1 . 1|) (including the 
chemical potential term) invariant. The minus signs which occur in ( |1.7|) are all 
important: We interpret all negative masses and densities as those pertaining to 
holes, or antiparticles. Thus, the duality transformations (|1.7| ) exchange particles 
and antiparticles. With this interpretation we always have 

^>0, (1.8) 
m 
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and thus the first two terms in ( |1.1|) are manifestly positive. (For more details see 
e.g. Section 3 of and references therein.) 

One of the important and interesting features of this collective-field formulation 
of the Calogero model is that it bears soliton solutions, which are known explicitly 
E' 12' 11' 1^1- They arise (to leading order in the expansion) as solutions of the 
classical equations of motion resulting from (|1.1|) [p!3l . 

Recently there has been considerable renewed interest in solitons in the Calogero- 
Sutherland models. In [|T^ it was shown that the collective-field theory of Calogero 
model is equivalent to a quantum version of an integrable Benjamin-Ono equation 



15]. In this way, it can be seen that the semiclassical one-soliton solutions studied 



m 



10| , |11|, |12| correspond to a single pole in a certain pole ansatz. 



The authors of [|^, have studied a specific duality-based generalization of 
the hermitian matrix model which is equivalent to a two-family Calogero model 
[p!8| , |T9|, The multi- vortex solutions of the coupled BPS equations were 

interpreted as giant gravitons [Q. In we have shown that the coupled BPS 
equations collapse effectively into a single-family BPS equation at special loci in 
parameter space of the two family model. A similar observation was also made 
concerning the more general coupled non-BPS variational equations. All this was 
the consequence of the invariance of the collective-field Hamiltonian of the two- 
family Calogero model under an Abelian group of strong-weak-coupling dualities 
[§], analogous to (pT?!) . In p4| and ||2^ it was shown that a large class of solitons 
in the two-family Calogero model can be obtained by reducing it to the effectively 
one-family Calogero model. 

These results on the collective-field solitons in various variants of the two-family 
Calogero model motivated us to revisit the collective-field theory of the original 
single- family Calogero model in [^, where we studied, among other things, the 
periodic soliton crystal solutions originally discovered in [^. In this paper we fo- 
cus on a particular type of such soliton crystals - the periodic BPS density waves 
of the collective field hamiltonian, and study their quantum stability. We present 
complete analysis of quadratic fluctuations around these BPS solutions. The cor- 
responding fluctuation hamiltonian (i.e., the stability operator) is diagonalized in 
terms of bosonic creation and annihilation operators which correspond to the com- 
plete orthogonal set of Bloch-Floquet eigenstates of a related periodic Schrodinger 
hamiltonian, which we derive explicitly. The resulting fluctuation spectrum is posi- 
tive, and therefore these density waves are stable. 

This paper is organized as follows: In Section 2 we review the solution of the 
static BPS equation associated with the collective Hamiltonian ( |1.1D . This is done 
by converting it into a Riccati equation which can then solve explicitly. The solution 
is a static periodic soliton - the finite amplitude wave solution of |10|. 

In Section 3 we study quadratic fluctuations around these BPS density waves 
and discuss their quantum-mechanical stability as was described in the paragraph 
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next to previous. The complete orthogonal set of Floquet-Bloch eigenstates of the 
related periodic Schrodinger hamiltonian, in terms of which we diagonalize the fluc- 
tuation hamiltonian, are derived explicitly in Section 4. Remarkably, the fluctuation 
spectrum is independent of the parameter which determines the density wave's am- 
plitude. As a consequence, the sum over zero-point energies of the field-theoretic 
fluctuation hamiltonian, and its ensuing normal-ordering and regularization, are the 
same as in the case of fluctuations around constant density background, namely, the 
ground state. We close by making some brief remarks concerning fluctuations around 
non-BPS density waves in Section 5. Finally, some technical details of the collective 
field formalism are relegated to the Appendix. 



2. Review of periodic BPS solutions of the collective field 
equation 

The Hamiltonian ( |1.1| ) is essentially the sum of two positive terms. Its zero-energy 
classical solutions are zero-momentum, and therefore time independent configura- 
tions of the collective field (pT3|) , which are also solutions of the BPS equation 

A - 1 ^ ^ ^ /■ dyp{y) 
P J x-y 

It is easy to check that the duality transformation (|1.71 ) maps a solution p{x) 
of (|2.1| ) with coupling A onto another solution p{x) = —Xp{x) of that equation with 
coupling A = As we shall see below in Eq. (|2.21|) , all solutions of are of 



B[p]^^^ + \4'^^ = 0. (2.1) 
2 p J 



definite sign, and never vanish along the real axis. Thus, such a positive solution of 



I]) is mapped by ( |1.7| ) onto a negative solution, and vice-versa. 



The BPS equation ( |2.1| ) may be written alternatively as 

(A - 1) d^p = 27rApp^ , (2.2) 
where p^ is the Hilbert-transform 

p«W = ir dy^ (2.3) 

of p. Note that for A = 1, where the CM describes non- interacting fermions, the only 
solution of ( f^.2D is p = po = const. This is also the case at the bosonic point A = 0. 
Henceforth, we shall assume A 7^ 0, 1. Space independent constant configurations 
p = Po are obviously solutions of (|2.1|) also for A 7^ 0, 1. However, for such values 
of A, ( |2.1|) bears also the periodic space-dependent density wave solutions, whose 
derivation we review in this section. 

The proper way to solve this nonlinear integro-differential equation is to consider 
it together with its Hilbert-transform[0, [1^, ^ 

(A-l)9.p^ = vrA((p^)2-p2 + p^), (2.4) 
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where on the RHS we used the identity''' 

2{pp^f = {py-p' + pl (2.5) 
(and the fact that dxP^ = {dxp)^ on the LHS). Here po is a real parameter such that 

dx (p(x) - Po) = . (2.6) 

It arises from the fact that we seek a solution of p(x) which need not necessarily decay 
at spatial infinity. Note that (|2.4|) is even in po. By definition, the sign of po coincides 
with that of p{x), the solution of ( p^ . A positive solution p{x) > corresponds to a 
BPS configuration of particles, and a negative one, to a configuration of antiparticles, 
as was mentioned following (|1.7|) . 

We proceed as follows. Given the density p{x), consider the resolvent 



oo 



— oo 

associated with it, in which z is a complex variable. 

The resolvent ^{z) is evidently analytic in the complex plane, save for a cut 
along the support of p(x) on the real axis. From the identity 

1 P 

±i7r5{x), (2.8) 



X ^iO X 
we obtain 

^±{x) = <^{x±iO) = p^{x)±ip{x). (2.9) 

Thus, if ^{z) is known, p{x) can be determined from the discontinuity of ^{z) across 
the real axis. 

An important property of ^{z), which follows directly from the definition ( |2.7|) , 

is 

oo 

QHz) = ^ [ fM^. (2.10) 

TT J \z-y\^ 

— oo 

Thus, if p(x) does not flip its sign throughout its support, we have 

sign (53 ^(z)) = sign (^3 2;)) sign (p(a;)) . (2-11) 

We shall use this property to impose certain further conditions on the solution of 
( Pjr^D below. 



t For a compendium of useful identities involving Hilbert-transforms see Appendix A of ||2^ and 
also Appendix A of the second paper cited in . 
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It follows from ( p. 91) that ( p.2|) and ( p.4|) are, respectively, the imaginary and real 
parts of the Riccati equation 



(A - l)9,$±(x) = 7rA($|(a;) + pi) (2.12) 

obeyed by both complex functions $±(x) . Let ^±{z) be the analytic continuations 
of $±(x) into the z— upper and lower half planes, respectively. These functions are 
evidently the two solutions of 

{\-l)dMz) = 7^K^{zf + pl). (2.13) 

subjected to the boundary conditions (^\{x+iQi) = $_(a;— iO) and sign + iO)) - 

sign {p{x)) = signpo, from (|2.9D. The resolvent ( ^.7] ) is then obtained by patching 
together ^+{z) in the upper half-pane and ^-{z) in the lower half-plane. 
The standard way to solve ( p.l3|) is to write it as 



^ ' d^<l>{z) = ik, (2.14) 



where 

*=^. (2.15) 

is a real parameter. 

Straightforward integration of (|2.14|) then yields the solutions 



2 _j_ ^ikz—u± 

M^)=^P0Y3-^, (2.16) 

where u± are integration constants. The boundary condition $^(a; -|- iO) = $_(x — 
iO) then tells us that m_ = —u'^. Clearly, '^u+ can be absorbed by a shift in x. 
Therefore, with no loss of generality we set 3m+ = 0. The second boundary condition 
sign (3$+(x + zO)) = signpo then tells us that u = 3?^+ > 0. Thus, ^±{z) are 
completely determined and we obtain (|2.7|) as 



_j_ „ikz— u sign (Qz) 

^^=^^0^3-^^—^. (2.17) 



As can be seen in (|2.21|) below, the density p{x) associated with ( p.l7| ) is indeed of 



definite sign, namely, signpo- 

The asymptotic behavior of ( 2.17| ) is such that 



$(±200) = ±ipo sign k . (2.18) 



This must be consistent with (|2.11 ), which implies (together with the fact that 



sign {p{x)) = sign po) that k must be positive. In other words, as can be seen from 
( p.l5|) , positive BPS density configurations (po > 0) exist only for A > 1, and negative 
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BPS densities (po < 0) arise only for < A < 1. The duality symmetry ( p..7|) , which 
interchanges the domains < A < 1 and A > 1, maps these two types of BPS 
configurations onto each other. Positivity of k, the condition ( p. .81) and the fact that 
sign p{x) = sign po imply that we always have 

signpo = signm = sign (A — 1) . (2-19) 

Now that we have determined ^{z), let us extract from it the BPS density p{x) 
and its Hilbert transform p^{x). From (|2.17|) we find that 



^ , , ^ , , — sin kx + i sinh u , , 

$+ X = $0r + zO = Po r ; , 2.20 

cosh u — cos kx 



from which we immediately read-off the solution of the BPS-equation as 

sinh u 

p{x) = Po 



cosh u — cos kx 



Tj . -. Sill IxiOC , . 

where both > and m > 0, and the sign of p(x) coincides with that of po. That p^ 
in ( p.21| ) is indeed the Hilbert-transform of p can be verified by explicit calculation. 
The static BPS-soliton, given by p(x) in ( p.21|) , is nothing but the finite-amplitude 



solution of 1 10]. It comprises a two-parameter family of spatially periodic solutions, 



all of which have zero energy density, by construction. The period is 

r = ^ = ^. (2.22) 

k Apo 

It can be checked by explicit calculation^ that 

^ j p{x)dx = po, (2.23) 

period 

oo 

and therefore that / (p(a^) — Po) dx = , as required by definition of po. Thus, 

— oo 

the parameter po determines both the period of the solution p(x), as well as its 
period-average, and the other (positive) parameter u determines the amplitude of 
oscillations between its extremal values 

u u 
Pmin = Po tanh - and pmax = Po coth - . (2.24) 



^Thc best way to do this computation is to change variables to t = e'*^^ and transform the 
integral into a contour integral around the unit circle. 
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Note also from (|2.23|) , that the number of particles per period is 



X 



X 



(2.25) 



A few limiting cases of (|2.21| ) are worth mentioning. Thus, if we let m — )■ , we 
obtain a comb of Dirac 5— functions 



p{x) 



A-1 



A 



6{x — nT) 



(2.26) 



If, in addition to m — )• 0, we also let k tend to zero (or equivalently, let the period T 
diverge), such that & = f remains finite, we obtain the BPS soliton solution pT| 



X 



X 



1 1 



6^ + 



(2.27) 



In fact, the original construction of the periodic soliton ( p.21| ) in |T0| was done by 
juxtaposing infinite solitons ( |2.27| ) in a periodic array. 

Note that the relation ( |2.25|) is preserved in both limiting cases discussed above, 
since the RHS of ( p.25| ) depends neither on u nor on k. 

Finally, by letting m — )■ oo in ( p.21|) , we obtain the uniform solution p = po of 
( p.l|) , independently of k. As can be seen from (|2.15| ), k blows up as A — )■ 1, namely, at 
the point of non-interacting fermions. Consequently, in this limit, the BPS density 
wave oscillates wildly between its extrema ( |2.24| ). This is clearly a pathological 
situation, unless these extrema coincide, which happens only when u — ?■ oo, namely, 
the uniform configuration p = po- The latter is the only solution of ( ^.l] ) at A = 1. 
Thus, if we seek match smoothly our rapidly oscillating solutions in the vicinity of 
A = 1 with the constant solution p = precisely at A = 1, we have to insist that 
the amplitude of the density wave should vanish as A — )■ 1 , i.e., let u — )■ oo. More 
quantitatively, according to ( p. 241) , a typical estimate of a derivative of the density 
wave ( p.21|) is k{pmax — Pmin) = J^hw; which for large k and u behaves asymptotically 
like fce~" ~ j^. Thus, we must have e" ^ in (|2.21| ), in order for the density 
waves to cross-over smoothly to the uniform solution at A = 1. None of these 
complications arise at A = 0, corresponding to non-interacting bosons. In this limit 



— )■ independently of u, and (|2.21|) tends to a uniform solution trivially. 



3. Fluctuations 

The collective-field formalism provides a systematic framework for the expansion 
^ By expanding the collective Hamiltonian ( |1 . 1|) around the static BPS so- 
lution in ( p.21|) , which we shall henceforth denote by Ps(x), we can go beyond the 
leading order and obtain the spectrum of low- lying excitations above Ps{x). Here we 
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shall concentrate on the next-to-leading terms, embodied in the quadratic fluctua- 
tions around Ps{x). 

To this end we write 

p{x) = psix) + r]{x) , (3.1) 

where ri{x) is a small density fluctuation around the wave solution ps{x), being 
typically of order Similarly, we shift p = ps + 6p, where /i^ = is the value of 

the chemical potential corresponding to Ps{x). Due to (|1.6| ), which is already satisfled 
by Ps{x), we clearly must have J dxrj = . This latter constraint is enforced by the 
shifted chemical potential 6p, as can be seen in ( |3.3| ) below. 

A convenient intermediate step is to expand the BPS combination B[p] around 
Ps{x). We obtain 



B[Ps + v] 



X 



Ps 



X 



ps 



+ ... 



(3.2) 



where the ellipsis stands for terms cubic in rj and higher, and where we used the BPS 
equation B[ps{x)] = 0, Then, substituting (|3^ ) in ([L1| ), and expanding Hcoii 



to second order in ri{x) , tt{x) and 6p, we obtain the quadratic fluctuation hamiltonian 
Hcoii^^^ around Ps{x) as 



H 



(2), BPS 



coll 



1 

2m 



2 1 

dxps{x) {dxTx{x)) + — 



dx Ps{x) 



2 \p. 



Xttt] 



H 



6p / dxT] + H, 



sing . 



(3.3) 



70) to eliminate all terms 



Here we used the constraint (|1.6|) and the BPS equation 
linear in the shifted quantities r] and 6p. 

Let us now evaluate -f^smg in the background of Ps{x). This will be also useful 



for later reference. The singular term multiplying p{x) in the flrst term in (|1.2|) is 
a constant, C. We can evaluate this constant as a divergent integral in momentum 



space. To this end, note that according to ( |3.15| ) below, C is proportional to the 
singular matrix element (a;||p||a;). Thus^ 



2m X — y 



Att 



y=x 



2m 



(X 



x) 



(A - l)k f dq 



4mpo 



(3.4) 



where we have also used (|2.15|) . Since C is constant, according to (p.6|) we may 
replace Ps{x) under the integral in ( p..2| ) by po- Therefore, Hging in this background 
is the same as in the familiar constant background p{x) = pQ. Combining J dxCpo, 



^This integral should, of course, be cut-off at some large momentum. The point is that the same 
regularization should be used for both the uniform and density wave backgrounds. 
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the first term in ( |1.2|) , with the momentum space representation of the second term 
there, we obtain our desired result as 



H. 



sing 



dx 



dq X-l 2 
271 4m 



k\q\) 



(3.5) 



We shall make use of this explicit form later. 

The hamiltonian ( p. 3D is essentially the sum of squares of local hermitian opera- 
tors, each of which is multiplied by the positive function Thus, the expectation 
value of ( p.3| ) with respect to any (normalizable) wave functional "^[rj] is positive, for 
all values of A. The fluctuation spectrum about Ps{x) is therefore positive, and the 
BPS density waves Ps{x) are stable. This, of course, comes at no surprise, since the 
energy density tied with Ps(x) is strictly zero - the lowest possible value for ( pTTI) . 

This positivity of (p. 3D, i.e., its quadratic structure, clearly calls for the intro- 
duction of the operator 



A{x) = dx7i{x) + i 
along with its hermitian adjoint 

A\x) = dx'rc{x) — i 



2 'Vp. 



Xnr] 



H 



(3.6) 



(3.7) 



It follows from the canonical commutation relations 



[n{x),T^iy)] = i5{x-y) , [ri{x),ri{y)] = [7r(x), 7r(?/)] = 
that A and A^ satisfy the commutation relation 

[A{x),A^y)]={l-X)d,dJ^-^^^]+2Xd, ^ 



Ps[X) 



x-y 



(3.9) 



with all other commutators vanishing. 

In terms of the operators ( |3.6| ) and ( p.7| ) , we may write the quadratic hamiltonian 
3f) as 

H^r' = ^ j dx ps{x) A\x) A{x) 

^ j P'^^^^ [A{x) , A\x)] + Hsing - 5p j dxTj. (3.10) 



As can be clearly seen by comparing ( p.9|) and (|1.2| ), the commutator term in ( p.lOp 
would cancel the singular term Hsing if Ps were constant, in accordance with 
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We shall see below (see (|3.21|) ) that this commutator term, in fact, exactly cancels 
the expression in (|375|) also in the background Ps{x). Thus, we obtain* 



H 



(2),BPS 
coll 



2m 



dx Ps{x) A\x) A{x) . 



(3.11) 



As we shall see below, this hamiltonian requires one last zero-point energy subtraction 
to render it finite. 

i'l) BPS 

3.1 Diagonalization of H^J^l in terms of bosonic creation and annihila- 
tion operators 

Our goal is to diagonalize the manifestly positive hamiltonian (|3.11|) . This could be 
achieved by decomposing -»/^^^A(x) and -v/^^^A*(x) into orthonormal modes. To 



this end, let us first observe that we may rewrite 



as 



where H is the single-particle hamiltonian 

1 



■ {x\H\ 



H 



:i-A) 



Ps{.x) 

2m 



P^^P-—\p\ 



Psjx) 
2m 



(3.12) 



(3.13) 



in which x and p are canonically conjugate position and momentum operators. The 
first term in ( p.l3| ), quadratic in p, can be read-off from ( |3.9[ ) in a straightforward 
manner. In fact, it has the standard non-relativistic hamiltonian form, since it is 
proportional to 



1 



-PVP^ 



P^ + ( ^d^\ogps{x) ] - ^d'^\ogps{x 



p^ + 



Ps[X) 

Po 



(3.14) 



where in the last step we used the BPS equation ( p.l|) . The second term in ( p.l3| ) 
arises because the matrix-element 



1 P 
{x\\p\\y) = -d^- 



which follows from the identity^ 



TT X — y 



\p\i){x) = -d^iijj^ix)) . 



(3.15) 



(3.16) 



■I- The term (5/i J dxrj merely constrains the zero- momentum Fourier mode of ri{x) to vanish, and 
we shall henceforth not write it explicitly 

^This identity can be easily established by applying the operator \p\ to the Fourier integral 
representation of tp{x), and using the formula (e*^^)-^ = le^^^signp. 
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In addition, in deriving the second term in (|3.13|) , one has to invoke equations (|2.1| ) 
and ( gjfSl) . 



The operator H is therefore a hermitian periodic Schrodinger operator, with 
period T = ^ (Eq.( p.2^ )). Consequently, it has a complete set of orthogonal Floquet- 
Bloch eigenstates ipq{x) and corresponding energy eigenvalues a;(g), where q is quasi- 
momentum. Note that H is not positive definite. Hence u:{q) may become negative 
over some range of q. 

As we shall show in Section 4 below, the spectrum of H has no gaps! This 
peculiar feature is obviously the result of the unconventional \p\ term in H . Since 
there are no gaps in the spectrum, we shall take the quasi-momentum q to range 
from — oo to +00 (i.e., work in the extended-Brillouine zone scheme). 

We shall now prove that the complete orthogonal set of Floquet-Bloch eigen- 
states ^q{x) can be used to diagonalize the fluctuation hamiltonian ( |3.11| ) over an 



appropriate Fock space. This can be done exclusively on basis of the general proper- 
ties of the periodic Schrodinger hamiltonian if, without any reference to the explicit 
form of the eigenstates ^q{x). We shall therefore defer derivation of the explicit form 
of the eigenstates ipq{x) and their corresponding energy eigenvalues uj{q) to Section 
4. 

Let us now list a few standard facts about the spectrum of H. By definition, the 
Floquet-Bloch eigenstates are quasi-periodic, and satisfy 

^q{x + T)=d'^^^q{x). (3.17) 

00 

The completeness relation / dq\uj{q)){uj{q)\ = 1, obeyed by the Floquet-Bloch 

—00 

eigenstates, may be written as 



dqipg{x)ipl{y) = S{x -y) . (3.18) 

Here, as usual, {x\u{q)) = (pq{x) . Orthogonality is expressed as 

00 

{uiq)\uiq')) = [ dx<f*(x)^gix) = 6iq-q'). (3.19) 



-00 



Finally, using the spectral decomposition H = f dqu{q) |ci;(g))(ci;(g)| we obtain the 

—00 

matrix element of H in the position basis as 

00 

{x\H\y)= [ dquiq)ipqix)ip*(y). (3.20) 
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Let us revert back to the discussion leading from (|3.1CI|) to (|3 . 1 1|) , namely can- 
cellation of (p.5|) by the commutator term in ( p. 10 ). By combining ( p. 12 ), ( 3.19|) and 
( p.20|) we can express that commutator term as 



/ dx Ps{x) \Aix) ^ A\x)\ = - [ dx (x\H\x) 

4m J ^ 2 J 



dxl / (3.21) 



In other words, this commutator term is nothing but the zero-point energy of a non- 
interacting bosonic field theory with eigenmodes given by the eigenvalues of H in 
( p.l3| ), w(g). We compute this dispersion relation in the next section. It is given by 
( [4.26| ). Substituting ( ^.26| ) in ( p.21| ) we see that it exactly cancels Hging in (p.5|), as 
promised. 

Since we seek a diagonalization of the fluctuation hamiltonian ( p.ll ) 



over an 



appropriate Fock space, we shall associate with each Floquet-Bloch eigenstate ipq{x) 
a pair of bosonic creation and annihilation operators a{q), a^{q) which obey the stan- 
dard bosonic commutator algebra 



[a{q),a\q')] = 5{q - q') , 
[a(g),a(g')] = [at(g) , at(g') ] = . 



(3.22) 



It is then a matter of straightforward calculation to show that the algebra (|3.9|) is 
realized by the normal-mode expansions 



y^A(x)= j dq\uj{q)\^^,{x) [e{u{q))a{q) + 9{-uj{q))a\q)\ 



-oo 
oo 



1 

y^^n^)= / dq\u{q)\^^l{x) [e{u:{q))a\q) + e{-uj{q))a{q)] 



(3.23) 



which take into account the fact that w(g) flips its sign. Indeed, it is easy to see that 
( ^.231 ) are consistent with [A(x), A(y)] = 0. In addition, fulfillment of ( p. 12] ), namely. 



2m 



2m 



{x\H\ 



(3.24) 



follows from (jpl), ( ^22]) , <^^, and (^^. Thus, ( ^231) is a legitimate 

realization of the operators A(x), A'^{x). Finally, upon substituting the normal mode 
expansions ( p.23|) in ( p.ll| ) we obtain, using the orthogonality relation ( p.l9|) , the 



- 13 - 



desired diagonal form of Hf^il^^^ as^ 

oo 

H^r' = I dq \u{q)\ [e (ooiq)) a\q)a{q) + 9 {-u{q)) a{q)a\q)\ . (3.25) 



This manifestly positive operator is not normal-ordered, due to the expansions ( |3.23| ), 
which contain both creation and annihilation operators. We should perhaps mention 
that a similar situation arises in the two dimensional anyonic model studied in [E7 



The divergent expectation value of ( p. 251) in the Fock vacuum is given by the sum 
over zero-point energies of negative u{q) modes, namely, 

oo 

(Oli^iir^O) = I |Kg)|^(-c.(g)) = -(|rfx) j |a;(g) 

— oo \q\>k 

= 2H,ing +ndx\ ^^k^ , (3.26) 



where in the last step we used ( |3.5| ) and ( |4.26| ). 

Note that the Floquet-Bloch eigenvalue u){q) given in ( [4.26| ) coincides with the 
dispersion relation of quasi-particles in the presence of uniform condensate po, namely, 
the ground state. Thus, according to (p.26|) , quantum corrections to the energy den- 



sity tied with the BPS-density wave studied here, due to quadratic fluctuations of 
the collective field, coincide with the analogous corrections around the uniform con- 
densate background. In other words, as one starts from the uniform solution of (|2.1|) , 
namely, ( p. 21 ) taken at infinite u, and then reduces u continuously to some positive 



finite value, there are no (next-to-leading order) quantum corrections to the energy 
density tied with the BPS-density wave relative to the ground state. (The classical 
values of energies in both cases, are of course null.) Moreover, since the single BPS 
soliton, or lump, is obtained as the limiting case ( |2.27| ) of the BPS-density wave, its 
mass it not corrected by quadratic fluctuations either. 

Finally, using the commutator algebra ( p.22|) to normal-order (|3.25|) , and sub- 



tracting the divergent contribution (|3.26|) , we obtain the desired diagonalized fluctu- 
ation hamiltonian simply as 

oo 

: {Hil]r') ■■= I dq \uj{q)\ a\q)a{q) . (3.27) 

— oo 

It is manifestly positive definite, as required. We need not worry about positivity of 
u{q) since only its absolute value enters ( ^.27| ). 



^Here we used 0'^{uj) ~ ^(w), as well as d{uj)d{—uj) = 0. The latter is responsible for the absence 
of aa and a^a^ terms on the RHS of (3.25). 



4. The exact Floquet-Bloch eigenfunctions and energy band 

Now that we have diagonahzed HfJl^^^ , it remains to determine the Floquet-Bloch 
eigenstates ^q{x) and corresponding eigenvalues u{q) explicitly, and establish their 
orthogonality and completeness. 

The eigenvalue equation H(j)q{x) = u{q)(j)g{x) may be written explicitly as 



where we have suppressed any g-dependence for brevity. Note also the slight change 
in notation of eigenstates from (fg to (pq. (This is done in order to avoid possible 
confusion in the discussion below, and we shall return to the original notation toward 
the end of this section.) This equation clearly calls for defining a new unknown 
function 

i'-^J (4.2) 
which has the same quasi periodicity ( p. 17 ) as (p. In terms of ip we have 

A\ 1 , „ , 



p^p_p^ + 27TXi^]\p\^|J = u^p. (4.3) 
2m J ps \2m/ 

Recalling ( p.l6|) and the BPS equation (p.l|) , we may write ( [4.3|) more explicitly as 
(A - l)dl^ - 2irX [pfd,^ + psd^t/j''] = 2moo^ . (4.4) 
4.1 Solving Schrodinger's Equation 

Equation ( |4.4| ) is very suggestive of taking its Hilbert transform, using the identity''" 

ifg"" + /^^7)^ = -f9 + fm , (4.5) 

where /o and 5^0 are subtraction constants of their corresponding functions, analogous 
to po in (|2.6|) . Since dx and the Hilbert transform commute, we find the Hilbert 
transform of ( [4.4|) as 



(A - 1)9^^^ - 27rA [pf 9,^^ - PsdxiJ + po^'o] = 2ma;^^ . (4.6) 



Life is made easier by the fact that ip^, the subtraction constant in (|4.6|), needed 
to render the integral of 'ip'{x) over the real axis, is null. This should be expected, due 
to Cauchy's theorem, under the assumption that ip{x) is analytic in at least one of 



^For derivation of (4^) for real functions see Appendix A of For a derivation of (4.5) 

for complex valued functions analytic in one half-plane, where they also decay at infinity (e.g., 
exponentials e*P^), see Appendix A in the second paper in [Q. If / and g are of the latter type, 
they are eigenfunctions of the Hilbert transform, with eigenvalue +i for analyticity in the UHP, 
and —i for the LHP. For such functions /o = go = due to Cauchy's theorem. 
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the half-planes where it also decays at infinity. We can also prove that ipQ = more 
directly. Since both iIj{x) and ip'{x) are quasi-periodic and obey similar relations like 
( ^TtD , it follows that 

°° ( ^ \ oo 



n=— oo 



= 27rV^(0)(e'«^-l)(5p(gT), (4.7) 

where 5p(x) is the periodic delta-function of period I-r. The last expression obviously 
vanishes, since '?/'(0) is finite by assumption. ^ Thus, t/'q = promised and it can 
be dropped from ( [4.61) . Reversing the argument, this independent proof that = 
implies that if tl}'[x) is analytic in one of the complex half-planes, it must decay to 
zero at infinity there. 

The natural thing to do now is to combine (|4.4|) and (|4.6|) , in a manner similar 



to (|2.12|) . Thus, adding and subtracting i times ( [4.4|) from (|4.6|) , we obtain the pair 
of equations 

(A - 1)9^F± - 27rA<l>±(9^F± = 2ma;F± , (4.8) 

where 

F^[x)=i)^[x)±ii){x) . (4.9) 
Recall that the two functions $±(a;), defined in ( p.9|) , are the boundary values 



$^(x) = $(2; ±iO) (4.10) 



of the meromorphic function $(-2;) in ( p.l7|) , as one approaches the real axis. Thus, 
$+(x) is analytic in the x— UHP, and <l>_(x) is analytic in the x— LHP. These are the 
coefficient functions in ( |4.8| ) . It thus follows from the theory of differential equations 
that -F+(x) should be analytic in the x— UHP, and that F^ix) should be analytic in 
the x-LHP. 

The two equations ( [4.8| ) can be united into a single equation 
(A - \)F'\z) - 2t:\^{z)F\z) = 2mujF{z) , 

or more compactly, 

F"{z)-^<!>{z)F'{z) = ^F{z), (4.11) 
Po A - 1 

for a single meromorphic function F{z), such that 

F±{x) = F{x±iO) . (4.12) 



•^Otherwise, if ?A(0) diverges, we can pick any other period of ipi^) s-t which endpoints ip is finite. 
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Our first step in solving ( [4 .111) is to remove the first-derivative term. Thus, 
following standard methods, we substitute 

F{z) = az)f{z) (4.13) 

where f{z) is a new unknown function, and ^{z) is determined by demanding that 
upon substituting ( [4. 13] ) in ( ^.llj ), the coefficient of f'{z) will be null. Thus, we find 

Using ( p.l7|) we may write this equation as 



(4.14) 



^1 j^J^ gjfcz—Msign (Qz) 
^ 2 ~^ 1 gjfc^— M sign (Srz) 



(4.15) 



The solution in each of the complex half-planes is immediate and is independent of 
the solution in the other half-plane. Thus, 



e = iC+ 9{Qz) 



ikz 

e 2 



-Qz) 



ikz 

e 2 



(4.16) 



with integration constants C±. In each half-plane, ^{z) is of course analytic, with 
discontinuity along the real axis. Granted more information on the analytic prop- 
erties of our solution ^{z) and F{z) we shall be able to determine these integration 
constants. 

We also need an expression for ^"/^ which appears in the equation for f{z). 
Taking the derivative of ( ^.14| ) and using ( |2.13| ) we find 

i! = ^ I^AV (2$2+ _ 

e 2po V 2po J \2poJ ^ ^ ^ 

With ( [4.17| ) at our disposal, we finally obtain the very simple equation 

2 



/"(--) + (5 

for f{z). We readily find the solution as 



^tpz 



(4.18) 



(4.19) 



with p a real momentum-like parameter.^ The resulting dispersion relation is, of 
course, 

A-1 



2m 



P 



(4.20) 



§More precisely, the two independent solutions of (i.l8) are e^'^^. Since they are related by 
flipping the sign of p, it is enough to consider only (4.19), since we allow p to range over positive 
and negative values. 
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Finally, combining ([4.19|) , (|4.16|) and ( [4 .131) we obtain the solution of ([4.11|) as 



iqz iqz 

F{z) = ei-sz) -— r— - 9{--sz) , (4.21) 



^ikz—u \^ J \ J ^ikz+u 



where 



\+P (4.22) 



is the quasi-momentum, as indeed, F{z + T) = e*''^F(z), in accordance with ( |3.17| ). 
In terms of g, the dispersion relation ( f4.2(J| ) may be written as 

^{q) = ^^(^ - q) ■ (4-23) 

Note that uj{q) is independent of the amplitude of Ps{x), which is governed by the 
parameter u. It depends on only through k, or equivalently, through po- 

We are not done yet, since we have to impose the condition that F±{x) = F(x± 
iO) be consistent with ( ^4.9| ) and also with the fact that V'o = ( [4.6| ), meaning that 
ip'{x) has to decay at infinity in the half-plane where it is analytic. Bearing in mind 
the last requirement, we obtain from ([4.9]) that 



F{x + iO) - F{x - 20) _ C+{q)e{q) e'"'' C_{q)e{k-q) e^-?^ 



2i 2 1 - e*''^-" 2 1 - 6*'=^+" 

,fr, , F(x + iO) + F(x-iO) tCJq)9(q) e^*^ iC(q)9(k-q) e^''^ 

ip [x) = = — 



2 2 1 ^ikx—u 2 \ ^ikx+u 

(4.24) 

It would be perhaps useful to spend a few words to clarify these equations. Let us 
concentrate on il){x). The first term in ip{x), analytic in the x— UHP, decays there 
at infinity only for g > 0, which is guaranteed by the 9{q) prefactor. Similarly, the 
second therm in ip{x), analytic in the x— LHP, decays there at infinity only for q < k, 
which is guaranteed by the other prefactor 6{k — q). With this assignment of domains 
of quasi-momentum, the first term in iIj{x) is an eigenvector of the Hilbert transform 
with eigenvalue i, and the second term is an eigenvector with eigenvalue —i, and 
therefore, ip^lx) in the second line of (|4.24| ), is indeed the Hilbert transform of tp{x). 



The domains of quasi-momentum in ( [4.24| ) appear asymmetric. In order to 
render them more symmetric, we shift q in the second terms in both lines of ( [4.24] ) 
by one unit of the reciprocal lattice, q = q' + k, where, of course, q' < 0. After a little 
algebra (including absorbing a factor — e~" in C_ and an appropriate redefinition 
thereof), and reinstating q for q', we can write iIj{x) more symmetrically as 

, cMOjq) e^'^ , C4q)ei-q) e^^^ 
'Vy^i - 2 1 - e^*:^-" 2 1 - e"**^^-" ' l^-^^J 
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With this new assignment of quasi-momenta, the dispersion relation ( [4 .231 ) may be 
written more symmetrically as 



^{q) = ^{k\q\-q'), (4.26) 



where g > corresponds to the first term in ( [4 .251 ), and g < to the second term. 
This function is made of two parabolas soldered together at g = 0. It vanishes at 
g = 0, ifc, with a cusp at g = 0. In addition, since > 0, it has two degenerate 
maxima at g = ±|, with maximal value Wmax = i^^^- Therefore, each value of uj{q) 
between w = and Wmax is four-fold degenerate, and occurs at < g < A; as well 
as at ±(/c — g) and — g . All other possible values of Ci;(±g), with |g| > /c, are just 
doubly-degenerate. 

We readily identify Co'(g) as the dispersion relation of quasi-particles in the pres- 
ence of uniform condensate po B HH ; as we have already mentioned following (|3.26|) . 



In particular, it is independent of the parameter m, which governs the amplitude of 
oscillations of the BPS-soliton. 

It is interesting to study the behavior of ( [4. 261) in the limit A — )■ 1. In this 
limit, according to ( p.l5| ), k diverges such that (A — l)k — > 27rpo, and therefore 
^{q) ^kl- This is, of course, the dispersion relation of sound waves propagating 
in the ground state of a one- dimensional Fermi gas of uniform density po and Fermi 
momentum kp = vrpo. This uniform Fermi gas configuration is the only solution of 
( p.l| ) at A = 1. Thus, the spectrum ( [4.26| ) of fluctuations around the BPS-sohton 
( p.21| ) behaves smoothly as A goes through A = 1. This supports our assertion, made 



at the end of Section 2, that these BPS-solitons should also be taken to be smooth 
around A = 1, which is achieved by setting e" ^ as A — > 1. 

4.2 The complete orthogonal set of Floquet-Bloch eigenfunctions 



Equipped with the explicit solution ( [4.25|) (and the definition (^4.2|) ), we can read-off 
the Floquet-Bloch wave-functions as 



A ( \ Po e'^- 

W) = \ — 7, TTl Ik^' 

" l-K Ps{x) 1 — e^"^ " 



27r Ps{x) I — e '■'^^ " 



such that 



— r-q ■ 



(4.28) 
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Here, the normalization factors C±. in ( [4 .251) are chosen such that (pq and cpq be 
normahzed to a delta-function, namely, 

dx (t)q{x)(j)*{x) = 5{q — p) q,p>0 



dx ^g{x)4>;{x) = 6{q - p) q,p<0, (4.29) 



which can be easily established by using the factorization 

1 pU 



Psix) 2posinhM 



Our work is not done yet, since (f)q and 0g do not form orthogonal sets. To see this 
choose some p > and g < and consider 

g~«(p+kl) I _ p-ikx-u 

2n 1 - 



M^mx) = , _ , (4.31) 



where we used ( [4.30| ) once more. The denominator can be expanded into a geometric 
series which we integrate term by term. Thus, we obtain 

/ dx (f)p{x)(j)*q{x) = J^e""" [S{p+ \q\ + nk) -e""5((n- l)k + p + |g|)] 

-oo 

= -e-"(5(p+ |g| - A;), (4.32) 

where in the last step we used the fact that k,p and |g| are all positive. In this region 
of parameters, the solutions ofp+|g| — A; = are any < p < and q = p — k (and 
therefore —k < g < 0). Thus, 0p(a;), with < p < A;, and (f)p-k{x) are not orthogonal, 
and we need to rotate them into mutually orthogonal combinations. This rotation 
will mix the two states, which is allowed physically, since both states are degenerate 
in energy according to ( |4.26D , and moreover, have the same quasi- momentum mod- 
k, i.e., both sates acquire phase e*^^ upon spatial shift by one period T. Obtaining 
these orthogonal combinations is a straightforward task, which essentially amounts 
to diagonalizing Pauli's matrix a^- The desired orthogonal combinations are found 
as 



,(+)/ N _ 'Pgjx) + 'i>q-k{x) _ psjx) , J{q-k)x\ 

^2(1 - e-«) " V 47rpo(l + e-) ^ 



V 2(l + e-) -V4vrpo(l-e-) " ^' 0<^<^- 

(4.33) 
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It is easy to check, using ( [4.32|) , that they are normahzed to a delta-function 

oo 

dx(f)^+\x)(f)'p+^*{x) = / dx(f)^f\x)(j)'f^*{x) = 6{q-p) 0<q,p<k. (4.34) 



By construction, they are also orthogonal to all the other eigenfunctions, with quasi- 
momentum larger than k or smaller than —k. In order to make our notations more 
symmetric with respect to assignment of quasi-momenta, let us rename in ( |4.33p 



'V(x) = -4-_;^ix) 



^iqx _ gi(g+fc)x^ ^ -A; < g < ^ 



(4.35) 

which of course, has no effect on the orthogonality relations. This concludes our 
derivation of the normalized orthogonal set of orthogonal Floquet-Bloch functions. 
For convenience, let us summarize them in the following list: 



4)q{x) 



2n ps{x) ' 



q > k 



Lfgix) = < 



i/2(l-e-") V 47rpo(l+e ") 



^2(l+e-") V 47rpo(l-e-") 



(^^iqx ^ ^iiq-k)x^ ^ <q<k 
PsM ^^iqa: _ g»(g+fe)x^ , -k < q < 



4)g{x) 



2n ps{x) l-e-'fc^-" ' 



q<-k, 



for which 

oo 

dx ipq{x)ip*p{x) = S{q-p) 
for all q,p E M, in accordance with ( |3.19D . 



(4.36) 



(4.37) 



4.2.1 The completeness relation 

We shall now verify that the functions in 
consider the LHS of (|3.18| ), namely. 



(^4.36| ) comprise a complete set. Thus, 



T{x,y) = I dq^l{x)^q{y) , 



(4.38) 
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and split it into contributions of the four sets of functions in ( [4.36D 



^i{x,y) = / dq(j)*{x)(f)q{y) 



dq4>\->{x) 4>[~\y) 



^i{x,y) 



dq(f)*(x) 4)g{y) 



(4.39) 



From the first and fourth hues in (|4.36| ) we obtain, by using the identities 



ikx 



X + ie 



p 

X 



-le 



ikx 



X — le 



that 



^i{x,y) 



l-e 



-2u 



iPol 



p 

-ie'^"" h it5(x) 

X 



1 



(4.40) 



j^^ik(y-x) 



2tx 



ps{x)ps{y) (1 



-ikx—u 



)(1 - e*'=2^-") y ~x + it 



T^{x,y) = Vl{x,y), (4.41) 
where, the last line follows from (|4.28| ). After some work we obtain 



Ti{x,y)+Ti{x,y) = 5{x-y)- 



^Ps{x)ps{y) sin^^ 



7r|po|(l-e-2") 5^ 



cos 



k{x - y) 



e "cos 



k{x + y) 
2 

(4.42) 



where we used the last equalities in (f4.40|) . In particular, the last cumbersome term 
in (|4.42| ) arises from the principal parts in ([4.41|) , which combine in such a way that 
the singular principal part is multiplied by a function which has a simple zero 
at X = y, which allows us to drop the P symbol. T2{x,y) and T2,{x,y) are evidently 
non-singular kernels. From the second and third lines in ([4.36|) , we can show, after 
some straightforward but tedious calculation that r2(x, y) + r3(x, y) exactly cancels 
the second, regular term in ( |4.42| ). Thus, 

T{x,y)=6{x-y), (4.43) 

proving completeness of the set of Floquet-Bloch functions ( |4.36| ). 
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4.3 Comments on uniqueness of the eigenstates 



Recall from the discussion following ([4.26|) that when —k < q < k, u){q) is four- fold 



degenerate. Thus, e.g., for < g < A;, also quasi- moment a k — q, —q and q — k all 
have a common value uj{q) (which ranges between w = and Wmax as Q varies in 
this domain). To this common eigenvalue correspond the four ort/io(7ona/ eigenstates 
(^^""(a;) and (i)^~\{x). This quartet of states splits into two pairs of 
states with common quasi periodicity, namely, (t)^(^\x) and which acquire 

phase e*"^^ under a shift of x by one period T (as can be seen from ( [4.36| ), in accordance 
with ( p.lTD ), and 0^^(x) and 0^^g(a;), which acquire phase e"*''"^. We can rotate each 
pair of these states by a 2 x 2 unitary matrix, which may even possibly be g-dependent, 
while leaving orthogonality (|4.37|) and completeness (|3.18| ) in tact. The rotated basis 



is as good as the original one given in i W^ . Thus, in this range of quasi-momenta, 
namely, —k < q < k, there is extra local energy dependent SU (2) x SU (2) symmetry. 
It would be interesting to investigate the origins of this symmetry further. 

In the remaining range of quasi-momenta, |g| > k, the energy spectrum is doubly 
degenerate, however, the corresponding two quasi-momenta ±g are not generically 
separated by an integer multiple of k, and these states cannot be mixed, except for 
a discrete set of values g„ = ±f A;, with n a positive integer, where there is an extra 
SU{2) symmetry. 

Finally, let us briefly comment on the zero-energy solutions. Consider approach- 
ing u = from within the domain of four-fold degeneracy. Let us pair the quartet of 
degenerate states according to their quasi-periodicities, as was discussed just above, 
namely, {(j)^f^\x) , {x)) on one hand, and (0q \x) , (j)l^\x)) on the other. These 
four limiting states are of course strictly periodic, and we can mix them by a uni- 
tary transformation. Thus, the energy dependent SU{2) x SU{2) is enhanced at 
w = to the much larger symmetry SU{4:). However, let us return to the paired 
states of the lower SU{2) x SU{2) symmetry. As can be seen from the explicit 
expressions in ( |4.36|) , we can take linear combinations of the members of each pair 
which are proportional to sj^^^- For such combinations, according to ( |4.2| ), we have 
ip{x) (X Ps{x). It is straightforward to check directly that ip^x) = Ps{x) is indeed a 
solution of ( ^73| ) when u = 0. Indeed, upon substituting these ip and u in ( |4.4| ) 
we obtain the derivative of the BPS equation ( p.l| ) multiplied by 2ps. One physical 
origin of this zero-mode has to do with the translational collective coordinate of the 
BPS density wave - i.e., it can be shifted arbitrarily in space^. However, due to the 
enhanced SU{4) symmetry, this cannot be the sole origin of zero eigenvalues in the 
spectrum of H, which requires further study. 



^Recall the imaginary parts 3u+ — which we absorbed as a shift of x following ( 2.16| ) 
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5. Conclusion and discussion 



In this paper we have completely diagonalized the hamiltonian of quadratic col- 
lective field fiuctuations in the background of BPS-density waves, which appear as 
static solutions of the equations of motion of the collective field formulation of the 
Calogero model. The fiuctuation spectrum is positive, demonstrating linear stability 
of the BPS-density waves. Remarkably, the fiuctuation spectrum a;(g) around these 
BPS-density waves coincides with that of fiuctuations around uniform condensates. 
The only difference between fiuctuations around these two background types is the 
explicit form of the orthogonal complete set of Floquet-Bloch mode functions (f)q{x). 
We computed these functions explicitly for the BPS-density wave background, by 
diagonalizing explicitly a related periodic Schrodinger operator, with a non-standard 
term, containing the absolute value \p\ of the momentum operator. Contrary to 
standard periodic Schrodinger operators, familiar from solid state physics, there are 
no gaps in the spectrum of the Schrodinger operator studied in this paper. 

We close this paper by making some brief comments concerning fiuctuations 
around non-BPS density wave solutions of the collective field equations of motion. 



These variational equations, for static configurations p{x), boil down top3 



B[p? - ^ d,{pB[p]) + 27rA - 2mp = , (5.1) 

P 

where the BPS combination B[p] was defined in Note that the BPS-density 

wave configuration Ps{x), given in ( p.21| ), is a solution of ( |5.1| ) corresponding to van- 
ishing chemical potential p = 0. In ||23[ we have found solutions of ( |5.1| ), which 
basically amount to subtracting from the BPS-density wave ( ^.21| ) either its min- 
imum or maximum value, thus obtaining extremal density waves which do vanish 
periodically, as opposed to the BPS-density wave. These solutions, which we named 
vortex crystals (corresponding to subtraction of the minimum, hence positive con- 
figurations ) and anti-vortex crystals (subtraction of the maximum, hence negative 



configurations), coincide with the large amplitude waves of in the static limit. 
We have discussed these vortex and anti-vortex crystals in |^3[ in detail, and com- 
puted the energy densities tied with them, which are all positive, hence above the 
zero-energy density tied with the BPS-density waves. Therefore the question of their 
linear stability poses an interesting problem, which is still open. 

This stability problem is more difficult to analyze than the stability of BPS- 
density waves, which we worked out here. To appreciate this difficulty, let us derive 
the quadratic fiuctuation hamiltonian in the background a non-BPS static solution 
p of (|5.1|). A convenient intermediate step is to expand 



B[p + r]] = B[p] + 



2 



(5.2) 



- 24 - 



where the elhpsis stands for terms cubic in r] and higher. Then, substituting ([57 
in and expanding Hcou to second order in ri{x) ,7i{x) and 6fi, we obtain the 

quadratic piece Hj:^^ as 



Hell = Vcou[p] - 5fl j dxr] + Hsing 

H / dx p{x) {dx7[{x))'^ -\ / dx p{x) 

2m J 2m J 



— [ dxB[p]r]T]", (5.3) 
m J 



where we used the constraint (|1.6| ) and the variational equation ( 5J. ) to ehminate all 



terms linear in the shifted quantities t] and 6fi. Here V^o««[p] is the energy tied in the 
extremal configuration p{x) (which was null in the BPS case). 

Now we can see the crux of the problem: p{x) is not a BPS configuration, and 
thus B[p] 7^ 0. This therefore generates a new term in (|5.3|) , which did not arise 
in the BPS case, namely, the last term in (|5.3|) , — (vrA/m) J dx B[p\ rjrj^ . This term 
prevents us from carrying out the factorization of the quadratic pieces in ( ^.31) in 
a manner analogous to ( p.ll| ), which proved so useful to the complete and explicit 
diagonalization of the quadratic fluctuation hamiltonian around BPS-density waves. 
A new idea is clearly needed to solve the non-BPS fluctuation spectrum. 



A. The collective field formulation of the Calogero model: 
derivation of Hcoii and its singular part Hging 

For the sake of being self-contained, we briefly review in this appendix the derivation 
of the collective field hamiltonian Hcoii in (pTll), and its singular part Hsing in ( |1-2|) . 



Standard references on the collective field formalism are [p!3| , p8| . Here we shall follow 
[ Pm , as well as @, which focus specifically on the Calogero model. 
The Calogero model, whose quantum hamiltonian is given by 

2m oxf 2m (xi — XjY 

i=l ' ij^j 

describes N identical particles of mass m living in one dimension. These particles are 
subjected to inverse- square pair interactions with dimensionless coupling A. ^ The 



^Note that we did not include in (AA) a confining potential. This is not really a problem, as 
we can always add a very shallow confining potential to regulate the problem (in the case of purely 
repulsive interactions), or else, consider the particles confined to a very large circle (i.e., consider 



(A.l) as the large radius limit of the Calogero-Sutherland model M). We shall henceforth tacitly 



assume that the system is thus properly regularized at large distances. 
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many-body wave functions are of the general form ^(xi, . . . , xjy) = A^S{xi, . . . , x^), 
where S{xi, . . . ,xj\r) is a function totally symmetric under any permutation of the 
particles, and A = nj<j(^* ~ -^i) Vandermonde determinant. The so-called 

Jastrow-factor A'^ arises due to the singular pair-interaction, which requires the wave 
function to vanish when any two particles coincide. The precise power-like vanishing 
is dictated by the requirement that the hamiltonian be self-adjoint. 

We see that for a generic value of A, under the interchange of any two particles, 
the wave function suffers a phase change of e*'^'^. These particles are therefore anyons. 
The singular pair- wise interactions vanish, of course, at A = 0, which corresponds 
to non-interacting bosons, and also at A = 1, which corresponds to non-interacting 
fermions. 

Powers of moments of the collective field p{x) in (|1.3|) are clearly the building 
blocks of all functions which are totally symmetric in the particle coordinates. Thus, 
the symmetric factor in the many-body wave-function, S{xi, . . . ,X]y), assuming it 
has a well-behaved large- limit, should become a well-behaved functional of p{x). 
It is precisely these symmetric wave-functions S{xi, . . . ,xn) on which the collective 
field operators, as well as the collective hamiltonian Hcou in ( [L.ID act. Thus, in order 
to transform ( [A. ID into a form amenable to collective field reformulation, we have to 



strip-off the Jastrow factors from the many-body wave function. This we achieve by 



performing on ([A.l|) the similarity transformation 

H ^ H = A-^HA^ , (A.2) 
It is straightforward to check that 

(E 5 - MA - 1) E - ^ E l-A-l- . V2 . (A,3) 

Thus, 

We can naturally interpret as part of a laplacian in some set of curvilinear coordi- 
nates which in addition to the Xi, also contain additional coordinates orthogonal 
to them. Vg is therefore the projection of the larger laplacian onto the subspace 
which depends exclusively on the original coordinates Xi. In particular, it is invari- 
ant under coordinate transformations which involve only the Xj. 

The space parametrized by the coordinates q"' is endowed with a metric Qab- We 
need not concern ourselves with the details of this metric (and of the additional 
coordinates, orthogonal to the Xj) except for the following two facts: First, we must 
clearly have-f 

V9 = A'\ (A.5) 
•l-The volume element must be positive. Thus, we should really interpret A^'^ — (A^)'^ ~ |A|^'^. 
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Second, in the subspace of the original coordinates obviously ds'^ = '^f^idx'^, and 
the corresponding block of Qat is therefore the unit matrix. 

At the particular values A = |, 1 and 2, these curvilinear coordinates are those 
of the symmetric spaces corresponding to real-symmetric, complex-hermitian and 
quaternionic-self-dual matrices, respectively, and is the projection of the laplacian 
into the singlet sector of the corresponding matrix space. 

By construction, H acts on the symmetric functions S{xi), and inner-products 
and matrix elements are computed with the integration measure 

dfi = A'^^d^x, (A.6) 

in respect to which H is hermitian. 

The collective field formalism amounts to performing a point canonical transfor- 
mation from the position operators Xi and their conjugate momenta pi = —i-^, 
to a new set of coordinates, namely, the collective field operators p{x) and their 
conjugate momenta n(x) = —ij^j^) and then expressing H in terms of p{x) and 
n(a;). 

A.l Point canonical transformations 

More precisely, we should look upon this as a special case of point canonical coordi- 
nate transformations 

q-^Q- = Q'^(q) (A.7) 

of the larger space, parametrized by the entire collection of coordinates g", in which 
only the subspace parametrized by the Xi is transformed into the new set of coordi- 
nates p(x), while the subspace orthogonal to the Xi remains unchanged. 
The metric in the new coordinates is given, in the usual manner, by 

and its inverse is given by 



Si.i,(Q) = 9.,.„(q(Q))i^i3^ (A.* 



SJ-(Q)=r"(q(Q))^g. (A.9) 

In these equations Q'"(Q) are, of course, the inverse coordinate transformations. It 



follows from ( [A. 81 ) that 

n{Q) = detnab = gJ^ (A.IO) 

where 

J = detf|^) (A.ll) 



is the jacobian of the transformation. Thus, we have. 



= y/gJ (A.12) 
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rendering the volume element invariant 



VUdQ = ,/gdci. (A. 13) 



Our next step is to transform the hamiltonian H in ( |A.4D to the new coordinates. 
Since Vg in (|A.3|) is that part of the invariant laplacian in our space, which in 
particular, remains invariant under transformations that change only the Xi while 
keeping the remaining orthogonal subspace unchanged, we simply have 



1 ^2 



2m 

where 



HQ = -—ns. (A.14) 



^ ^Qab^A. ) , (A.15) 



^ VndQ-\ dQ\ 

and the subscript s in Vq^, means a projection on the subspace originally parametrized 
by the Xj. By construction, of course, Vg^ = V^. 

This new expression for H = Hq is symmetric with respect to the measure 
a/cJc/Q. Life would be much easier if we could rid ourselves of this potentially com- 
plicated measure, and map H onto an effective hamiltonian -ffg// which is symmetric 
with respect to the fiat measure (iQ. This we can achieve by performing the similarity 
transformation 

H,ff = Q^HqsQ-'^ . (A.16) 

Let us now massage Heff into a more transparent form. In order to avoid 
cluttering of our formulas, we shall henceforth not display the subscript s explicitly. 
Thus, all the formulas below should be properly projected onto the subspace originally 
parametrized by the Xj. It is a matter of straightforward calculation to show that 

^ 9g» V dQ^ J 

^ r^'^'^Jn- - (^^-{^^'0,)) - C^n^'C, (A.17) 



where we have defined 



dQ^ \dQ'' 



Ca = ^(\0gQ),a (A. 18) 



and where (■)^a indicates a derivative with respect to The operator f24 Vq f2 4 
is manifestly symmetric with respect to the fiat measure (iQ, as is evident in each of 



the lines in (|A.17|) . Thus, Heff is indeed the desired hamiltonian we set out to find. 



which, following ([A.17|) , we may write explicitly as 



1 h"^ 

Heff = — Pa n"'' Pb + — Ca n"'' Cb + — (n"'' Cb) (A. 19) 

2m 2m 2m ^ ^ ' 
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where we introduced the momentum operators 



P, = -ttg^. (A,20) 



and displayed h dependence exphcitly. The terms in ( A.19| ) quadratic in h may 



be thought of as a generahzation of the centrifugal barrier which arises in the ra- 
dial hamiltonian in D dimensions^. Evidently, these terms are purely a quantum 
mechanical effect. 

It is easy to see that 

Ca = (A.21) 

where F^^ is the second Christoffel symbol (i.e., the connection) of Qat- However, 
sometimes a direct computation of the Ca from their definition (|A.18|) , or from the 



identity ( |A.21| ), may be too difficult to carry in practice. Thus, in order to bypass 
these potential difficulties, we shall now derive an identity satisfied by the Ca, from 
which we could compute them with somewhat less effort. 

To this end we argue as follows: The invariant laplacian acting on a function 
which is a scalar under coordinate transformation produces yet another scalar func- 
tion. Thus, 

VjV'(g) = V|^(g(Q)) . (A.22) 

In particular, the coordinate functions themselves are scalars (their differentials are 
one-forms). Let us define the quantities 

= - WjQ" = -hVlQ" . (A.23) 



It follows from the definitions ( |A.23| ) and ( |A. 18| ) that 



or 



+ 2n n^'Cb + h b = , (A.24) 



which is the desired identity to determine the Ca- To simplify the computation, 
we are free to choose in ( |A.23| ) the coordinates in which the computation of 
u"" = —KV'^Q'^ is as simple as possible. In our case, these are just the original 
coordinates. 



§The radial part of the D-dimensional laplacian — r^^^^^^ dr{r^^^ dr), defined with re- 
spect to the measure r^^^dr, may be transformed by a similarity transformation into Wf. ~ 
r"^2^ V^r~"^2^ = dl — '^^^^l^.F^'^^ = df. — Cr^^^Cr, which is defined with respect to the fiat 
measure dr. 
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A. 2 Application of the point canonical transformations to the collective 
field problem 

The transformation from the Xj's to the density field p{x) makes sense only in the 
limit — 7- oo, since on one side there are A^ position operators Xj, and a continuum of 
density operators p{x) on the other side. The continuum of these density operators 
are not all independent. For example, they are subjected to the constraints ( |1.5| ) 
and ( |1.6| ). Thus, in order to facilitate this transformation, one has to regularize the 
continuum theory. This is most conveniently achieved in momentum space. As our 
independent collective variables we choose the first A^ Fourier modes 



N 

p^= dxe-'^''p{x) = e-'^""^ (A.25) 

i=i 



cut-off at some kmax, where k is properly discretized, e.g., by putting our system in 
a large box of size L, and imposing appropriate boundary conditions. The details of 
this discretization are not important for our discussion of the large- A^ limit. Assume 
next that the particles condense and that the mean particle density in the condensate 
is of the order of some value p. Thus, the microscopic inter-particle distance will be 
of the order ^ ~ ^. Consequently, the maximal Fourier components should be of the 
order kmax = y ~ P- Thus, the high density limit makes k^ax (and letting 

L — 7- oo in the end, makes k continuous). The collective field reformulation of the 
model is therefore valid in the high density regime, where the system behaves like a 
continuous medium. 

For large but finite spatial box size L, momenta are discrete. We then take the 
first lowest A^ Fourier modes pk as our new coordinates Q"^. ^ Since the metric gab 
in the original coordinates is given by ds'^ = YliLii^x^Y + . . ., i.e., its block in the 
subspace of interest is simply the unit matrix, we obtain from ( |A.8D that 



fi''-^*' = n{k, k'; [p]) =J2—^ = -kk'pu+w (A.26) 

1=1 ' * 

in obvious notations, where momentum modes play the role of the new coordinate 
indices. 

The effective hamiltonian ([A.19|) contains also the quantities Ca, which we will 



determine from the identity ( |A.24D . Thus, we have to compute u^'' = u{k; [p]). From 
the definition (|A.23|) we obtain" 

u;(A;;[0]) = -V^0fc, (A.27) 



was stressed above, this transformation affects only tlie subspace originally parametrized by 
the Xi. The coordinates of the orthogonal complement subspace remain unchanged. 
"From this point on we set H — 1 again. 
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where was defined in ( |A.3|) . Thus, we obtain (see Eqs.(3.11) -(3.13) in 0) 



Adding and subtracting a Yli term in the second term in ( |A.2S ) leads us to the final 
expression 



u{k; 



(1 - \)k'^pk + AA;^ / dapkapk{i~a) ■ 



(A.29) 



We should now substitute ( [A.26| ) and ( |A.29D in ( |A.24| ). It is easy to see that the 
last term there vanishes: 



(A.30) 



Thus, we obtain from ( [A.24|) that 

u[k; [p]) + 2Y,Q{k,-k'- [p])C{k';[p])=0 



(A.31) 



which we can use to determine C[k; [p]). In the combined limits of large density and 
infinite spatial box, the fc-sums tend to Fourier integrals. Thus, in the limit, 

/rik' 
— nik,-k';[p])Cik';[p])=0, (A.32) 

which we shall now transform to x— space. To this end we need (see Eq.(6) in ) 
as well as 



f e^'Mk; [p]) = (A - l)dlpix) + 2XdJ ^^^M^ dy , (A.33) 
2tc J x-y 



dkdk' 
(2^ 



^ikx+ik'y 



n{k, k'; [p]) = d.^dy {p{x)6{x - y)) 



(A.34) 



Using ( |A.33D and (|A.34|) , we Fourier transform ( |A.32| ) to x-space and obtain 
uj{x; [p]) - 2d,^ {p{x)d^C{x; [p])) = 

2d^ 

2 J X — y 

which we can readily solve for d^C, and obtain 



^ ^d,p{x) + X-[ ^^^^ dy - p{x)d,C{x; [p]) 



(A.35) 



dxC{x; 



A - 1 d^p 
2 p 



+ A 



pjy) 

x-y 



dy. 



(A.36) 
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We now have all the ingredients required for computing the effective hamiltonian 
( [Al9l) . Thus, substituting ( [OiD , ( [OeD and Pa := n(x) = in the first two 



terms in ( [A.19|) we readily obtain the first two terms in ( |1.1|) . The singular piece 
Hsing in ( [1.2D arises from the last, divergence term in ( [A.19|) . From ( |A.24| ) and from 
( |A.30|) we thus have 



-L(r^-a) =-,^-\:=-,^ (A.37) 
2m ^ 4m 4m J 5p{x) ^ ' 



from which if^mg in (|1.2|) follows. 
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